† Corresponding author. E-mail:
Project supported by the Russian Foundation for Basic Research (RFBR) (Grant Nos. 16-32-00043 and 14-02-97026).
The effects in electrostatic models of chevron surface-stabilized ferroelectric liquid crystals are investigated through numerical modeling. To study smectic C* director distribution within the cell, we consider two nonlinear approaches: the chevron interface does not interplay with the electric field; the electric field interplays with the chevron interface. The obtained results of the director field distribution are compared with the earlier linearized studies. We find that whether or not the electric field interplays with the chevron interface, the electro–optic response requires a generalized approach for its description. The threshold electric field, which is necessary for switching between two stable director states in the chevron cell is evaluated. This study suggests that, in many cases of practical interest, electro–optic response to the electric field and the threshold electric field can be precisely estimated. We argue that, beside being numerically efficient, our approach provides a convenient and a novel standpoint for looking at the electro–optic response problem.
Considerable efforts have been made since Meyer et al’s pioneering work[1] to exploit the physics of ferroelectric liquid crystals (FLCs) in smectic C* (SmC*) phase. Correspondingly, a number of studies in SmC* liquid crystals confined between parallel boundaries and subjected to homogeneous boundary conditions is now available. For the so-called “bookshelf” geometry, the orientational order of the long SmC* molecular axis in thin films exhibits a one-dimensional positional order of the molecules arranged in a layered structure.
In the “bookshelf” geometry, smectic layers are normal to the cell surfaces, and often spontaneously deform into the chevron structures due to strong positional anchoring. In such structures, the layer edges at the boundary surfaces are not shifted, but the layers are tilted, and meet in the center of SmC* sample, which is called the chevron tip. Consequently, the understanding of the chevron structure is of considerable importance for the development of applications based on FLCs.
A series of theoretical models of chevron structures has been developed in recent decades. The pioneering theoretical description of the chevron structures in terms of a discontinuity in the layer tilt, while maintaining a continuous director structure at the chevron interface was given by Clark and Rieker.[2] A model of compressible layers was considered in studies.[3–5] Another approach assumed that the molecular tilt in FLCs is related to the layer tilt.[6] In Ref. [7], Limat extended the model,[4] where unequal cone and layer tilt angles are allowed. A similar phenomenological Landau–de Gennes approach was employed in a series of other studies (see e.g., Refs. [8]–[11]).
The early studies[2,12] consider models, where the layer tilt angle is fixed, i.e., ±δ0 in each part of the chevron, but the azimuthal angle φ is varied between the plates, while the cone angle θ is the constant. Thus, the system represents two flat crossing planes at the angle 2δ0. Continuous researchers’ interest to FLCs generated another model with fair simplifications,[13] where Nakagawa’s approach[3,4] was developed.
The chevron structure of smectic layers is a well-known characteristic of most surface-stabilized cells filled with FLC in the SmC* phase. These cells generally exhibit two stable director states. They can be switched between the stable states by an external voltage applied to the cell. Switching occurs only at fields stronger than the threshold field Eth because the director rotation is accompanied by a significant increase in the elastic energy around the chevron tip.[2] The bistability of such cells does not allow getting the desirable gray scale performance.[14,15]
Field-induced molecular reorientations in a monolayer of SmC* can be investigated by determining the spatial dependence between the azimuthal angle φ of the director and the applied electric field. Earlier, simplified models of biaxial modeling[16] and electro–optic response[13,17,18] were developed. The main simplifications in these studies regard the assumption that the layer tilt angle and the cone angle are small, which correspond to the temperature range close to the SmC*-SmA* (smectic A*) phase transition that has many dislocation defects.[19] Therefore, for lower temperatures[20] the electro–optic response of SmC* is expected to be different from the previous models with fair simplifications, but for a narrow temperature range.
In this study we aim to develop the generalized electrostatic model of the chevron FLC in the SmC* phase. On the basis of the continuum theory, we derive the functional, which represents the total free energy density of SmC*, and explains why the dielectric anisotropy term for weak electric fields should not be ignored. To calculate the spatial dependence of the azimuthal angle, we minimize the free energy density assuming that the smectic layers interplay with the electric field and do not. The threshold field Eth is calculated for the generalized model.
The plan of the article is as follows. The model and basic equations are presented in Section 2. In Section 3 we present numerical simulations under the assumption that the chevron configuration does not depend on the electric field. In Section 4 we turn our attention to the model when the chevron interface couples with the electric field. Details on computation of the threshold electric field are sketched in Section 5. Concluding remarks are given in Section 6.
We examine a cell of thickness d filled with SmC* in the chevron configuration, and confined between two parallel substrates. Let the cell be maintained at a temperature far below the SmC*–SmA* phase transition.[19,20] It is known that molecules of SmC* form a layered structure and are tilted from the layer normal
The total free energy density f of the layer includes the elastic free energy density felas, contribution to the electric energy density felec, and deformation of smectic layers fu, i.e., f = felas + felec + fu. A commonly used term for felas in the one-constant approximation can be written as
The spontaneous polarization vector
In order to introduce the contribution to the free energy density, associated with mechanical dislocations of smectic layers, consider a molecular length a0 and the smectic layer spacing h along the z axis. The layer spacing can be found by the identity
With the above assumptions, the total free energy density takes the form
Supposing that the shape of the chevron layer displacement field u(x) does not depend on the electric field, we minimize f with respect to φ. After routine trigonometric simplifications, the Euler–Lagrange equation takes the form
The chevron interface in the cell center impacts on the orientational ordering of the director field
To substantiate the assessment of the φ-dependencies with common parameters d = 5 μm, θ = 22.5°, K = 5 × 10−12 N, εa = 2.5, P0 = 2.9 × 10−4 C·m−2, we plot Figs.
As we expected, consideration of common polar angle and space-dependent layer tilt angle yields several differences in spatial distribution of the c-director field from its profile based on the linearized models for small δ and θ. Our computations indicate that in the generalized approach the φ-dependence differs from that found from models.[13,18] The proposed model indicates that weak electric field is unable to reorient the c-director parallel to the plates due to the interplay between the P0E term, strong boundary conditions, and non-coupling layer tilt with the electric field (Fig.
When the director is in the U-state, the deviations between the generalized model and the simplified model are less significant for the weak electric field and the strong electric field because the c-director has already oriented almost parallel to the boundary plates (as shown in Fig.
A number of factors can influence the layer configuration when SmC* molecules interact with the electric field. At this point it is worth examining the possible layer and director configurations within the cell. The system will show equilibrium structures of the c-director and layer configurations that reach minimum states of the total free energy density (
In Fig.
Layer bending, which is characterized by δ(x) is greatly dependent on the constant B.[13] As indicated in Fig.
To find the threshold field Eth, we suppose that the director switches from the D- to the U-state, i.e., rotating by angle π. When the director switches to the U-state, molecules rotate along the allowed trajectory forming the structure, which is energetically unfavorable, and determines the potential barrier W [the smectic cone angle reduces to zero, i.e., the liquid crystal in the smectic A phase]. Therefore, the difference between free energy of the system in the state Fπ/2 with φ(0) = π/2 and free energy of the D- or U-state F0 represents the potential barrier, i.e., W = Fπ/2 − F0.
Having determined profiles for δ and φ from the model represented in Sections 3, the electrostatic forces holding SmC* dipoles at their equilibrium state can be expressed by the integral
The results obtained in this study are more general because we studied non-linearized potential (
We have presented a theoretical study that determines the spacial dependence of the director field by taking into account a detailed layer structure and common parameters of SmC*. Many recognized electro–optical phenomena in SmC* with a chevron structure could be satisfactorily explained by this generalized description. Hereby, the presented study of electro–optic response in FLCs can be successful when we cannot ignore a number of terms in Eq. (
Two separate aspects of the influence of the electric field on the layer structure were studied. First, the model of non-coupled smectic layers with electric field. Second, the problem of coupled layers with the electric field. The latter problem represents the system of nonlinear differential equations of the second order, and it gives comprehensive results. Calculations of the azimuthal angle distribution in such layers were investigated by numerical methods within the framework of the continuum model. The obtained δ- and φ-dependencies were compared with the results derived from the model for small layer tilt and cone angle. The differences between the models of our focus imply that the issues of the study may be useful in optimization of mesomorphic properties of liquid crystal cell.
The threshold electric field, which is necessary for switching between two stable configurations was determined for non-coupled electric field with smectic layers. The generalized model is consistent with the previous studies, and confirms that Eth is inversely proportional to the cell thickness. To calculate Eth for the model of coupled electric field with smectic layers, a more accurate description of the interaction between the chevron structures and the electric field is required.
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 | |
11 | |
12 | |
13 | |
14 | |
15 | |
16 | |
17 | |
18 | |
19 | |
20 | |
21 | |
22 | |
23 | |
24 | |
25 | |
26 | |
27 | |
28 | |
29 | |
30 | |
31 | |
32 |